Microcavity polariton-like dispersion doublet in resonant Bragg gratings 
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Periodic structures resonantly coupled to excitonic media allow the existence of extra intragap 
modes ('Braggoritons'), due to the coupling between Bragg photon modes and 3D bulk excitons. 
This induces unique and unexplored dispersive features, which can be tailored by properly designing 
the photonic bandgap around the exciton resonance. We report that ID Braggoritons realized 
with semiconductor gratings have the ability to mimic the dispersion of quantum-well microcavity 
polaritons. This will allow the observation of new nonlinear phenomena, such as slow-light-enhanced 
nonlinear propagation and an efficient parametric scattering at two 'magic frequencies'. 



Since the pioneering proposals by Yablonovitch and 
John [T], photonic crystals (PhCrs), structures charac- 
terized by a spatially periodic dielectric function, have 
attracted enormous attention due to their rich physics. In 
particular, the occurrence of photonic bandgaps (PBGs), 
i.e. frequency regions where the propagation of light is 
strongly inhibited and the extraordinary ability to ma- 
nipulate and control the photonic flow, make PhCrs very 
appealing for many applications py. Furthermore, the 
current accurate engineering of photonic states allows 
to investigate light-matter phenomena in PhCrs; one of 
the most fascinating being the strong-coupling regime 
of light-matter interaction [5] . Non-trivial modifications 
of the photonic dispersion are expected to occur when 
PBG materials are coupled to polarizable excitonic me- 
dia PlH], as in this case the true eigenmodes of the mixed 
photonic-excitonic system are exciton-polaritons [5], i.e. 
the normal modes born from the strong coupling between 
Wannier-Mott (WM) excitons [3, 4 and Bragg photons 
propagating in the PBG structure. 

In this Letter, we show through theoretical analysis 
that, by embedding a three-dimensional (3D) exciton res- 
onance within the PBG of a one-dimensional (ID) Bragg 



grating, the dispersion relation of the resulting exciton- 
polariton states can, under certain conditions, mimic the 
dispersion of a doublet of microcavity (MC) polaritons, 
i.e. the quasiparticles resulting from the coupling be- 
tween two-dimensional (2D) excitons in a quantum well 
and the optically confined photons in a semiconductor 
planar microcavity [5 . The very peculiar dispersion fea- 
tures of these Bragg polariton modes have not been un- 
veiled previously, and, as we demonstrate here, give rise 
to a wide variety of new nonlinear phenomena, absent 
in standard gratings, such as slow- light-enhanced nonlin- 
ear propagation and an ultra-efficient parametric scat- 
tering at two 'magic frequencies'. In addition, owing to 
the extremely small effective mass of these novel hybrid 
exciton-photon modes, new routes for the appearance of 
macroscopic coherence in solid state systems are opened. 

It is possible to rigorously show, starting from 
Maxwell's equations, that the interaction between the 
electromagnetic field in a periodic medium (like the 
ID grating here examined) and the exciton-polarization 
wave, is described (in dimensionless variables) by the fol- 
lowing polaritonic coupled-mode equations (PCMEs): 



l[dr+%h]f+ld<:f + Kb+Pf/2 = 0, (1) 

t[dr + %i,]b-id^b+nf+pf,/2^0, (2) 

idrPf + (hx + d)pf + [\pf\^ + 2|pbP] Pf + Pf^ 0, (3) 

idrPb + (i7x + d)pb + [\pb\^ + 2\pf\^] Pb + pb = 0, (4) 



In Eqs. (T]|4l, the fields {/, b} represent the amplitude of 
the slowly- varying envelopes of the forward and backward 
propagating electric field respectively, while {pf,pb} are 
the corresponding quantities for the exciton polarization 
field. Spatial dispersion terms have been omitted for sim- 
plicity, as their contribution to the propagation dynamics 
is negligible. The dimensionless temporal and longitudi- 



nal spatial variables are r = t/to and C = z/zq respec- 
tively, with to = Cb /wq and zq = eta I \f^^ ^b being the av- 
erage background dielectric function of the grating, and 
ujQ is the central pulse frequency. The ID grating is de- 
scribed by the dielectric function e(z) = eb{\ + ^[e*'^^^ + 
g-ifcBzjy'2|^ where /i ^ 1 is the depth of the grating, 
and kB = 2J^loqIc is the grating Bragg wavenumber. 
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The dimensionless frequency detuning between the Bragg 
frequency lob and the exciton resonant frequency Wx is 
given by d = AutQ, where Aw — ojb ~ i-^x- The hn- 
ear absorption of the background medium 7ph and exci- 
ton oscillator damping 7x (i.e. the exciton homogeneous 
linewidth) are also considered in our model through the 
dimensionless quantities 7ph,x = io7ph,x- The polariton 
splitting LOc, which measures the interaction between the 
transverse components of the excitonic polarization and 
the retarded electromagnetic field, is specified by the di- 
mensionless parameter p = tQUj'^ / (2ujj^) . Finally, solution 
of the linearized Eqs. (l]|4l yields the dispersion rela- 
tion of Bragg polaritons 6 = d{q), where 6 = Aw/wq 
and q = Ak/ko are the frequency and wavenumber de- 
tunings from the bandgap center respectively, and have 
been introduced via a phase shift of the fields of the kind 
exp(iqC — iSr). 

Throughout the whole paper we use a design based on 
Zinc Oxide (ZnO) as a representative example of our the- 
oretical calculations. ZnO has received substantial atten- 
tion in recent years as a material for blue light emitting 
devices, transparent electrodes, and solar cells. This is 
due to the robustness of its exciton resonances, which 
remain stable up to room temperature, and its excep- 
tionally large binding energy and oscillator strength |13j . 
Following recent advances in ZnO growth and deposi- 
tion, clean exciton resonances have been produced using 
molecular beam epitaxy [T?], pulsed-laser deposition [1], 
sputtering [TB] and led to the realization of a ZnO MC 
Hj. We use ZnO as our excitonic material [see also Fig. 
(sFe)], and we therefore design the photonic crystal in or- 
der to have the PBG located around uj^. The second kind 
of layer is made of Zr02. The free FX^^^(r5)-exciton 
binding energy of ZnO is hujc « 60 meV, and the exciton 
central frequency is hio,. « 3.3771 eV, which corresponds 
to Ax — 367.4 nm [2]. The linear absorption coefficient 
and the exciton homogeneous linewidth of ZnO depend 
considerably on the fabrication technique. However, at 
T = 5 K we make the estimates 7ph ~ 2.5 meV and 
7x ~ 0.25 meV [T5]. The strength of the exciton-photon 
coupling is thus p/to = uj1/{2u}^) ~ 5.33 x 10""* eV. Near 
the exciton transition frequency, ZnO has a refractive in- 
dex Til ~ 2.37, while for Zr02 n2 — 2.326 at the same 
wavelength. Thus we have the following grating parame- 
ters, to be used in Eqs. (l]|4l: p = An^/eb ~ 0.0187, and 
K = peb/4: ~ 0.0258 ^ 1, satisfying the shallow grating 
condition on which our PCMEs are based. The quarter 
wavelength condition gives the layers' width: Li — 38.8 
nm and L2 — 39.9 nm for ZnO and Zr02 respectively. 
The nonlinear refractive index of ZnO at Ax is estimated 
to be nNL(ZnO) « 10"" mVW at 5 °K, about 9 orders 
of magnitude larger than that of bulk silica [5] . 

Equations (l]|4l represent the foundational result of 
this Letter. Previous attempts to describe polaritonic 
gratings have been made in the framework of the linear 
transfer matrix method (TMM) [3 , approaches that do 
not take into account the large nonlinear optical response 
of excitons near resonance, or by using a model based on 
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FIG. 1: (Color online) (a,b,c) Polariton dispersion on the 
[5, g)-plane for three different values of d. Dashed lines 
are the bare Bragg-photon and exciton dispersions. Blue 
solid lines indicate upper-photon (UP) and lower-photon (LP) 
branches, red solid lines are upper-Braggoriton (UB) and 
lower-Braggoriton (LB) branches, (d) Reflectivity Rj (blue 
solid line) and absorbance Aj (black dotted line) as functions 
of 5 for a finite grating (L ~ 14 pm for our proposed ZnO 
design). The dispersion is shown with red dashed lines. 



Maxwell-Bloch equations [7^ , but without taking into ac- 
count the crucial importance of both forward and back- 
ward propagating exciton polarization waves. Our ap- 
proach provides a simple, complete and tractable set of 
equations that avoids the complications of the TMM, but 
at the same time yields results that are very close to the 
ones obtained by numerically solving the full Maxwell 
equations coupled to the material equations. 

Figs, [ija-c) show the evolution of the coupled and 
bare linear dispersions when varying d from vanishing 
to negative values. Coherent coupling between 3D bulk 
excitons and Bragg photons gives rise to two new intra- 
gap modes (Bragg polaritons, or Braggoritons [3J), the 
dispersion of which can be efficiently tailored by simply 
tuning d, for instance by slightly modifying the grating 
parameters around the fixed exciton resonant frequency. 

The four branches in question [which we label upper- 
and lower-photon (UP,LP) and upper- and lower- 
Braggoriton (UB,LB) branches, see Figs, [ija-c)] have 
a striking resemblance with the dispersions of a doublet 
of coupled MC polaritons [8]. As shown in Fig. [ijc) for 
a detuning d = —0.04, the pair (UP,LB) form the first 
half of the doublet (named MCI), while its counterpart 
(named MC2), formed by the pair (LP,UB), is reversed, 
and not placed symmetrically with respect to MCI for 
d^O. 

Eqs. (l]4) allows calculation of the linear refiectiv- 



ity Rj and absorbance Aj as a function of d for finite 
gratings, see Fig. [ijd). Physical dispersion branches are 
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FIG. 2: (Color online) (a) Physical parts of the dispersion 
branches (for which R{5) < 1) are indicated with blue solid 
lines, for d — —0.04. Red dots indicate the position of the 
two inflection points on the UB and LB branches, which are 
located at q's that have the same magnitude but opposite 
signs, see also Fig. jsjb). (b) GVD [d'^qj /dS^) as a function 
of (5. Red dots indicate zero-GVD points. (c,d) SoHd blue hnes 
are the group velocities Vgj/c — {dqj /d5)~^ as functions of 
q, for d = and d = —0.04 respectively. Red dashed lines 
indicate the dispersion. 
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FIG. 3: (Color online) (a) Effective polariton mass parameter 
Mj as a. function of d. (b) Schematics of the parametric am- 
plification process occurring in the proximity of the inflection 
points of both Braggoriton branches, for d — —0.04. Grey 
areas indicate the photonic sub-gaps, (c) Matrix element for 
the intraband parametric scatterings shown in (b), for the LB 
branch (blue line) and for the UB branch (red line), (d) Same 
as (c) but for the LP branch (blue line) and for the UP branch 
(red line), (e) Sketch of our proposed design for a polaritonic 
Bragg stack. The excitonic material is ZnO (shown in red. 
Ax = 367.4 nm), while the second material is Zr02 (shown in 
blue). 



selected with the criterion R < 1, see also Fig. IJa). It IS 
evident that the polariton feature introduces an anoma- 
lous transmittance region located inside the bandgap, at 
a detuning d from the bandgap center, which divides 
the bandgap into two smaller sub-bandgaps [3 . Fur- 
thermore, we observe here that the formation of such 
sub-gaps is accompanied by the appearance of two zero 
group -velocity dispersion (GVD) points, therefore divid- 
ing the spectrum into alternating regions of anomalous 
[{d'^qj/dS'^) < 0] and normal [{d^qj/d6^) > 0] GVD [see 
red dots in Figs. [2ja,b)]. It is worth noticing that in 
conventional gratings such points do not exist, while MC 
polaritons possess only one inflection point on their lower 
branch, which is of paramount importance for the occur- 
rence of the nonlinear effects observed in experiments 
[HI [H] ■ We shall see that the full significance of this anal- 
ogy goes well beyond the mere similarity in the dispersive 
properties. 

Due to its large excitonic component and flat disper- 
sion, the UB (LB) branch with d < (d > 0) turns out 
to be a potential candidate for slow-light-enhanced non- 
linear optics at low intensities [10'. This is confirmed by 
calculating the group velocities (GVs) of all branches for 
d = and d = -0.04 [Fig. ^c) and (d) respectively]. 
The UP/LP branches show GVs approaching the speed 
of light in the medium [normalized to ±1 in Fig. [2]^c,d)] 
for large values of 6, as in conventional gratings. The 
GV can be reduced considerably when 5 approaches the 
band edges, but then progressively less photons will be 
available due to evanescence of the electric field in those 
regions [10 . However, Figs. E[c,d) show that in proxim- 
ity of the UB /LB branches (0 ^ d) , small GVs can be 
obtained. For d = 0, Fig. [2|c), the central GV curve is 
symmetric, while, for d ^ 0, the GV curve is asymmet- 
ric, see Fig. [2j^d). Close to d ~ d, exciton absorption can 
strongly affect both UB/LB branches, depending on the 
precise value of 7x. Being mostly excitonic in nature and 
spectrally narrower, the UB (LB) branch is the one most 
affected by exciton absorption for d < (d > 0). More- 
over, a reduction in group velocity comes at the price of 
a reduced bandwidth, which is a fundamental limitation 
for all slow-light devices [10 . This corresponds here to 
the progressive straightening of the fiat UB branch for 
large values of |d|. At low temperatures, typically T < 5 
K, there are frequency regions for which absorption is 
reasonably small and relatively far from the band edges, 
where the group velocity is greatly reduced. For instance. 
Fig. [2][d) shows that an increase of the group index up 
to 20 is achieved for realistic parameters corresponding 
to the ZnO-based grating structure discussed above. By 
increasing the interaction time between the medium and 
light field of approximately an order of magnitude, it is 
possible to achieve low- light- level nonlinear optics |10j . 
This will open up new vistas for slow-light-enhanced non- 
linearities in gratings. 

We have calculated that the physical effective polariton 
masses at g = for each dispersive branch mj are given 
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by 



(5) 



where Mj{d) = [{d^Sj /dq'^)]~lQ is a d-dependent coeffi- 
cient shown in Fig. [sj^a). The prefactor fnr ^/et, / {2cX^) 
assumes the value of w 4 • 10~^me for bulk ZnO, which 
is of the order of magnitude of the cavity photon mass 
in the upper polariton branch of MCs, the latter being 
in practice not usable for nonhnear optics purposes due 
to the absence of inflection points in the polariton dis- 
persion [8]. Thus ID polaritonic gratings allow the pos- 
sibility to obtain very small polariton masses compared 
to the lower branch of conventional MC-polaritons (the 
mass of which is ~ lO^^We ^). The above consider- 
ations, valid for ID structures, suggest that when the 
concept of resonant polaritonic grating is extended to a 
2D periodic system (by considering for example struc- 
tures similar to those described in Ref. [H]), there will 
be realistic chances to obtain dispersions possessing ex- 
tremely small masses, thus suitable for experiments on 
Bose-Einstein condensation (BEC) and superfluidity of 
polaritons [12j. 

Owing to their excitonic content, Braggoritons are very 
nonlinear due to the strong interactions of their con- 
stituents, thus giving rise to strong optical nonlinearities 
of Kerr (cubic) type. The nonlinear coefficient at the ex- 
citon resonance can be several orders of magnitude larger 
than the one for highly-nonlinear optical fibers [6]. We re- 
port here that a stimulated parametric scattering process 
occurs for Braggoriton states, which is the analogue of 
the parametric amplification in semiconductor MCs origi- 
nally proposed and demonstrated in 2000 by Savvidis and 
Baumberg In this experiment, when an intense pump 
excites the lower MC polariton branch at a certain 'magic 
angle' (corresponding to a magic wavenumber, located 
near the inflection point), one observes a large amplifl- 
cation of a weak probe beam which stimulates the same 
MC lower branch at normal incidence. Such an amplifica- 
tion is due to the scattering of the two pumped polariton 
states into a pair of signal and idler polaritons [9], which 
follows from the phase-matching conditions (energy and 
momentum conservation). Only modulationally unsta- 
ble frequencies can grow rapidly enough to produce the 
exponential amplification of the weak probe (located at 
q = 0) at the expense of the pump. In our system, the 
phase-matching condition is provided by the nonlinear 
terms in Eqs. ([3]|4]), and is satisfied in close proximity of 
two inflection points [located in the UB and LB branches 
as shown in Fig. |2j^a,b)] - note that only one inflection 
point is available in MC polaritons - thus effectively giv- 
ing rise to amplification of Braggoriton states at g = 0. 
As we are dealing with a grating structure, these two in- 
flection points can be externally excited just by changing 



the frequencies of the input pump and probe pulses, see 
also Fig. |3jb) and (e). Hence, Braggoriton amplification 
is characterized by two magic frequencies [5mi.m2 in Fig. 
[sjb)], in correspondence of which one transfers polaritons 
from the pump {q — q„i) to states with q — Q [see Fig. 
|3jb)]. The final states are located symmetrically in both 
frequency and wavenumber space. 

The scattering processes on the various dispersion 
branches can have very different matrix elements, and 
thus the efficiency can vary considerably for different 
branches, as shown in Figs. |3jc,d). The intraband ma- 
trix elements of the UP and LP branches are shown in 
Fig. [3](d). They are negUgible for small values of \d\, 
as they are not phase-matched, and thus the efficiency 
of the process on these branches is very low. The op- 
posite scenario is observed for the UB/LB branches, see 
Fig. |3jc). As a typical example, when d ~ —0.04, at a 
value of \q\ k, 0.38, where the zero GVD points are lo- 
cated, the ratio between the matrix element for the UB 
branch and the one for the LB branch is approximately 
i^UB/LB = \Mub/Mlb\ ~ 6, so that the UB branch is six 
times 'more nonlinear' than the LB branch for the spe- 
cific material parameters chosen. This ratio can be tuned 
up by several orders of magnitude by just increasing the 
value of \d\. For instance, for d = —0.2 we would have 
i^UB/LB ~ 2 X 10'^ at the zero GVD points. This allows 
a very large tunability of the effective nonlinearity of the 
various branches in polaritonic gratings. 

In conclusion, we have studied some previously unex- 
plored nonlinear properties of ID Bragg gratings, when 
the PBG is near-resonant with an excitonic feature of 
the medium. We have developed a new model based on 
a set of coupled- mode equations [Eqs. (l]|4l] that is able 
to accurately describe the linear and nonlinear propa- 
gation of Braggoriton states, which will be instrumental 
for future theoretical investigations. The novel feature of 
such gratings is the formation of two inflection points in 
the dispersion characteristics, leading to the formation 
of regions of slow-light enhanced nonlinear propagation, 
as well as the existence of a strong parametric ampli- 
fication at two magic frequencies. We believe that the 
results shown in this Letter are particularly important in 
the field of nonlinear parametric processes with exciton- 
polaritons. We are confident that the present work will 
stimulate further investigations of similar effects in PhCr 
structures of reduced dimensionality (for both exciton 
and photon fields) with potential applications to polari- 
ton quantum optics and polariton macroscopic coherence 
in mixed excitonic-photonic systems. 
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